Entanglement increase from local interactions with not-completely-positive maps 
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Simple examples are constructed that show the entanglement of two qubits being both increased and decreased 
by interactions on just one of them. One of the two qubits interacts with a third qubit, a control, that is never 
entangled or correlated with either of the two entangled qubits and is never entangled, but becomes correlated, 
with the system of those two qubits. The two entangled qubits do not interact, but their state can change from 
maximally entangled to separable or from separable to maximally entangled. Similar changes for the two qubits 
are made with a swap operation between one of the qubits and a control; then there are compensating changes of 
entanglement that involve the control. When the entanglement increases, the map that describes the change of 
the state of the two entangled qubits is not completely positive. Combination of two independent interactions that 
individually give exponential decay of the entanglement can cause the entanglement to not decay exponentially 
but, instead, go to zero at a finite time. 
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I. INTRODUCTION 

We construct simple examples here that show the entan- 
glement of two qubits being both increased and decreased by 
interactions on just one of them. In our first and basic step, 
taken in Sec. II, we have one of the two qubits interact with a 
third qubit, a control, that is never entangled or correlated with 
either of the two entangled qubits and is never entangled, but 
becomes correlated, with the system of those two qubits. In 
Sec. Ill, we do this for each of the two entangled qubits, and 
consider the combination of the two interactions, with sepa- 
rate control qubits that are not correlated and do not interact 
with each other The two entangled qubits do not interact, but 
their state can change from maximally entangled to separable 
or from separable to maximally entangled. Similar changes 
for the two qubits are made with a swap operation between 
one of the qubits and a control; then there are compensating 
changes of entanglement that involve the control. This is de- 
scribed in Sec. II. A. 

Whenever the entanglement increases, and in some cases 
where the entanglement decreases, the map that describes the 
change of the state of the two entangled qubits is not com- 
pletely positive and does not apply to all states of two qubits. 
It all depends on whether there are correlations with the con- 
trols at the beginning of the interval for which the dynamics is 
considered. The maps are described in Sec. IV and discussed 
in Sec. V. The completely positive maps that decrease the en- 
tanglement have akeady been described 1 1]. 

When the interaction of each qubit with its control by itself 
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gives exponential decay of the entanglement, the combination 
of the two interactions gives exponential decay at the rate that 
is the sum of the rates for the individual interactions, when the 
two interactions are made the same way. Making them differ- 
ently can cause the entanglement to not decay at that rate or 
at any single rate. Instead, the entanglement goes to zero at a 
finite time; the state becomes separable and remains separable 
at later times. This is described in Sec. III. A. Similar behavior 
has been observed in more physically interesting and mathe- 
matically complicated models |2, 3, 4]. 

These examples are built on the same framework, but to a 
very different design, from those we made for Lorentz trans- 
formations that entangle spins |5]. There the momenta that 
played the roles of controls were purposely coiTelated. Here 
the controls are kept independent. The framework makes the 
operations transparent by describing the qubit states with den- 
sity matrices written in terms of Pauli matrices, so you can 
see the Pauli matrices being rotated by the interactions. States 
are shown to be separable by writing out the density matrices 
explicitly as sums of products for pure states. For each in- 
teraction here, the map that makes the change of the density 
matrix for the entangled qubits is described by a simple rule 
that particular Pauli matrices in the density matrix are multi- 
plied by a number; equivalently, the map of the state of the 
entangled qubits is described by a rule that particular mean 
values are multiplied by a number. 

Our examples show that statements like "entanglement 
should not increase under local operations and classical com- 
munication" 01 are not generally true outside the set of 
local operations considered in the original proof In our 
examples, each control qubit interacts with only one of the 
two entangled qubits. In this sense, the quantum operations 
are local. CoiTelation with a control at the beginning of the 
interval for which the dynamics is considered can give local 
operations that increase entanglement. 
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II. ONE INTERACTION 

We consider the entanglement of two qubits, A and B. We 
use Pauli matrices Si, E2, S3 for qubit A, and Pauli matrices 
Si, S2, S3 for qubit B. We let qubit A interact with a third 
qubit, which we call C. We think of C as a control. By inter- 
acting with qubit A, it will control the entanglement of qubits 
A and B. We work with states represented by orthonormal 
vectors |a) and \(3) for C. We consider a state of the three 
qubits represented by a density matrix 



n 



1 



■if 



(2.1) 



For p± this gives 

P± = Q ± (SiCos</> + S]2sin0)Si 



1 \ 1 



2/4 



± (-Si sin + S2 cos 0)^2 - S3S3] 
[1 ± (Si cos (/i — S2 sin 0)Si 

± (Si sin + S2 cos <j>)E2 - S3S3] 



- -[]I±(SiSi+S2S2)c0S<^-S3S3] 

= p±cos^{(j)/2) + p^sin2(0/2). 



(2.8) 



with p the density matrix for the state of qubits A and B. 

We follow common physics practice and write a product of 
operators for separate systems, for example a product of Pauli 
matrices S and S for qubits A and B, simply as SS, not Scg)S. 
Occasionally we insert a (g) for emphasis or clarity. We write 
^A, 1b, Ic^ but we do not put labels A and B on the Sj and 
Sfc. The single statement that the Sj are for qubit A and the 
Sfc are for qubit B eliminates the need for continual use of 
both A and B lalels and signs. 

Suppose p is one of the density matrices 



P± = ^(1 ± S^iSi ± S2S2 - S3S3). 



(2.2) 



Both p+ and p_ represent maximally entangled pure states 
for the two qubits. They are Bell states. The state of zero total 
spin is represented by p_ and the state obtained from that by 
rotating one of the spins by tt around the z axis is represented 
by p+. 

For a rotation W, let Da{W) be the 2x2 unitary rotation 
matrix made from the S, so that 



Da{W)^ ^Da{W) = iy(S) 



(2.3) 



where M^(S) is simply the vector S rotated by W. Let W{4)) 
be the rotation by <j) around the z axis, and let Da{4i) be 
DA{W{cj))). 

We consider an interaction between qubits A and C de- 
scribed by the unitary transformation 



U ^ DAma){a\ + DAi-(f>mm 



or, in Hamiltonian form. 



U 



-i<j>H 



(2.4) 



(2.5) 



with 



i7 = S3-(|a)(a|-|/3)(/?|). (2.6) 

This changes the density matrix p for qubits A and B to 

p' - i:rc[{U®lB)Tl{U ®Ib)^] 

= ]^Da{(1>)pDa{c^)'' + \da{^<I>)pDa{-c^)^. (2.7) 



A. From maximally entangled to separable and back 

We focus first on the case where is 7r/2. Then both p+ 
and p_ are changed to 

P' - J[1-S3S3] 

^(1-1^3)^(1+ S3) 

^) ^(1+^3)^(1-^3). (2.9) 

The density matrix p for a maximally entangled state is 
changed to the density matrix p' for a separable state that is 
a mixture of just two products of pure states. The inverse of 
the unitary dynamics of qubits A and R takes p' back to p; it 
changes a separable state to a maximally entangled state. 

The dynamics continuing forward also changes this separa- 
ble state to a maximally entangled state. As goes from 7r/2 
to TT, the density matrix p'j_ changes from that of Eq. ( 12.91 ) to 



P± = Pt- 



(2.10) 



There can be revivals of entanglement between two qubits 
when there is no interaction between them, as well as when 
[8] there is. 

Changes in the state of qubits A and B from maximally 
entangled to separable and back to maximally entangled can 
also be made very simply with a swap of states 1 9] between A 
and C. This can be done with a unitary operator U with 
U a unitary operator for qubits A and C that acts on a basis 
of product state vectors simply by interchanging the states of 
A and C. There is interaction between qubits A and C only; 
qubit B is not involved. 

Applied to an initial state described by Eqs. ( 12.1b and ( 12.21) . 
where qubits A and B are maximally entangled, this swap op- 
eration gives a separable state for A and B. Applied a second 
time, it restores the initial state where A and B are maximally 
entangled. For qubits A and B, this is similar to what hap- 
pens when goes from to 7r/2 to tt. For the three qubits, 
it is different. The swap operation does not change the com- 
plete inventory of entanglements for the three qubits. It just 
moves the entanglements around. In particular, C becomes 
maximally entangled with B. We will see, in Sees. II. C and 
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D, that the interaction described by Eqs. ( I2.4l i. i2.5\ and i2.6\ 
does change the complete inventory of entanglements for the 
three qubits. When the state of qubits A and B changes from 
maximally entangled to separable and back to maximally en- 
tangled, there are no compensating changes of other two-part 
entanglements. In particular, qubit C never becomes entan- 
gled with anything. 



B. Concurrence 



The change of entanglement is smaller when (f> does not 
change by tt/2. ^From Eq. ( I2.8l l we have 

= i[l±(SiSi + S2S2)cos0-K(I]iSi)(S2S2)], (2.11) 

after rewriting the last term. This shows that for both p'_^ and 
p'_ the eigenvalues are 



), -(1-COS0), 0, 



(2.12) 



because SiSi and S2S2 each have eigenvalues 1 and —1 and 
together they make a complete set of commuting operators: 
their four different pairs of eigenvalues label a basis of eigen- 
vectors for the space of states for the two qubits. The Wooters 
concurrence jlOtl is a measure of the entanglement in a state 
of two qubits. It is defined by 



Cip) 



0: 



Ai — v A2 — v A3 — vA 



(2.13) 



where p is the density matrix that represents the state and 
Ai,A2,A3,A4 are the eigenvalues, in decreasing order, of 
p S2S2 p* S2S2, with p* the complex conjugate that is ob- 
tained by changing S2 and S2 to — E2 and — S2. From 
Eq. (12.11b we have 

p'± S2S2 {p'^r S2S2 = p'± {p'±r {^2^2? ={p'±f (2.14) 

so for the \fXl are the eigenvalues of and the concur- 
rence is 



(2.15) 



We can consider the change of entanglement as changes 
through any interval. When | cos 0| decreases, the entangle- 
ment decreases. When | cos 01 increases, the entanglement 
increases. 



C. Two-part entanglements 

The only two-part entanglements are when qubit A is in 
one part and qubit B is in the other There is entanglement 
between qubit A and the subsystem of two qubits B and C and 
between qubit B and the subsystem of two qubits A and C, as 
well as between qubits A and B. There is never entanglement 



between the state of qubit C and the state of the subsystem of 
two qubits A and B. The density matrix 

{U®lBmU®tB)^ = \DA{4>)p±DA{4>)^\a){a\ 

+ lDA{~c^)p±DA{-4>)^mm 

(2.16) 

is always a mixture of two products of pure states. The re- 
duced density matrix for the subsystem of qubits A and C, 
obtained by taking the trace over the states of qubit B, is just 
Ia ® 1 c/4, and the reduced density matrix for qubits B and 
C, obtained by taking the trace over the states of qubit A, 
is ts ® lc/4. There is never entanglement or correlation 
between qubits A and C or between qubits B and C. The 
reduced density matrices for the individual single qubits are 
just 1^/2, 11b/2, and lc/2. The only subsystem density ma- 
trix that carries any information is the density matrix p for 
the qubits A and B, which is changed by the interaction with 
qubit C. The entropy of the subsystem of qubits A and B can 
increase or decrease, but there is no change of entropy for any 
other subsystem or for the entire system of three qubits. 



D. Three-part entanglement 

There is three-part entanglement. The state represented by 
the density matrix (12.16b is called biseparable because it is 
separable as the state of a system of two parts, with C one 
part and the subsystem of two qubits A and B the other part. 
It is not separable as the state of a system of three parts A, B, 
and C. The density matrix ( 12.161 ) is not a mixture of products 
of density matrices for pure states of the individual qubits A, 
B, and C. If it were, its partial trace over the states of C, the 
reduced density matrix that represents the state of the subsys- 
tem of the two qubits A and B, would be a mixture of products 
for pure states of A and B. That happens only when cos (f>isO. 
In that case, we can see that the density matrix ( 12.161 ) is not a 
mixture of products for pure states of the individual qubits A, 
B, and C because its partial transpose obtained by changing 
S2 to — S2 is not a positive matrix. 

In the classification of three-part entanglement for qubits, 
biseparable states are between separable states and states that 
involve W or GHZ entanglement fUl [H [H. Let Hi, 
112, be Pauli matrices for the qubit C such that |a)(a| 
is (1/2)(1 +03) and \f3)(f3\ is (1/2)(1 -03). Bounds from 
Mermin witness operators say that for separable or bisepara- 
ble states 



-2< (E,s,n,- 



-SjSfcnfc-EfeSjnfc-SfeSfen,) < 2 (2.17) 



for j,k — 1,2,3 and j ^ k; a mean value outside these 
bounds is a mark of W or GHZ entanglement ITill . In our 
examples, these mean values are always 0. A mean value 
{\GH Z){GH Z\) larger than 3/4 for the projection operator 
onto the GHZ state. 



\GHZ) 



1 



V2 



|0)|0)|0) 



1 



i>|i)|i), 



(2.18) 
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is a mark of GHZ entanglement; it can not be larger than 3/4 
foraiy state [1 3.1 . A mean value {\GHZ){GHZ\) lai-gerthan 
1/2 is a mark of a VF state; it can not be larger than 1/2 for 
a biseparable state (vi\. In our examples, {\GH Z){GH Z\) 
is always 0. A mean value {\W) (W\) larger than 2/3 for the 
projection operator onto the W state, 



\W) 



1 



V3 



1>|0>|0) 



ID 



1 



1), (2.19) 



with a rotation around the x axis. If neither cos (1)a nor cos 4>b 
is 1, the concurrence becomes zero, the state separable, be- 
fore cos(/)^ or cos (t>B is zero. The interactions of qubits A 
and B with their controls change maximally entangled states 
to separable states. The inverses change separable states to 
maximally entangled states. In the following subsection, we 
describe the density matrices that show explicitly that the sep- 
arable states are mixtures of products of pure states. 



is a mark of W entanglement; it can not be larger than 2/3 for 
a biseparable state li 1 3i] . In our examples. 



1 



{\W){W\) = -(l±cos 



(2.20) 



This mean value does not involve either entanglement or cor- 
relation of the qubit C; it would be the same if both |a) (a| and 
1/3) (/3| in the density matrix ( 12.16b were replaced by {l/2)c, 
the completely mixed density matrix for C. 

For any 0, the density matrices ( |2.16l l for the two cases + 
and — are changed into each other by the local unitary trans- 
formation that changes the Pauli matrices for one of the qubits 
A or i? by rotating its spin by tt around the z axis. As a func- 
tion of 0, the mean value (|Vl^)(Vl^|) changes in opposite di- 
rections for the + and — cases. So will any mean value for the 
states described by the density matrices ( 12.161 ). if it changes at 
all. 

For the states described by the density matrices ( 12.161 ), the 
only nonzero mean values that involve the qubit C are 



(SsSiHa) = ±sin( 



(2.21) 



These would be the same if they were calculated with only the 
|a) (q;| part or only the |/3) (/3| part of the density matrix ( 12.161 ). 
In fact, they are the same as (SiS2)(n3) and (I]22i)(n3) 
calculated for one of those parts. Their values do not requke 
either entanglement or correlation of C. 



III. TWO INTERACTIONS 

If a control were coupled similarly to qubit B as well, then 
cos (j) would be replaced by cos (pA cos (/i^ in the next to last 
Une of Eq. ( 12.81 ) and in Eqs.( |2.11| ) and ( 12.151 ). If the coupling 
of qubit B is made with a rotation around the x axis instead 
of the z axis, then the next to last line of Eq. ( |2.8l l becomes 

p± = -[1 ± SiSi COS0A ± COS0^ COS0B 

-SsSscosc^b]. (3.1) 

Rewriting the last term and looking at eigenvalues in terms of 
SiSi and S222 as before yields the concurrence 

C(p^)=^ max[0, I cos0yi| + | cos cos + | cos 0_b| — 1]. 

(3.2) 

When cos 4>a is 1, these Eqs. dS.ll l and ( 13.21 ) describe the result 
obtained when there is only the interaction of qubit B made 



A. Exponential decay 



To describe exponential decay of entanglement we let 



cos0^ = e COS0B 



(3.3) 



by letting each interaction be modulated by a time-dependent 
Hamiltonian H{t) that is related to the Hamiltonian H of 
Eqs.dlS) and by 



H[t) = H-^ = HTcotcf), 



(3.4) 



where (j) and F are (pA and F^ or and F b ■ The same result 
could be produced in different ways. The interactions could be 
with large reservoirs instead of qubit controls 12. SH]. Each 
qubit ^ or i? could interact with a stream of reservoir qubits 
11511. Here we are interested in the way the entanglement is 
changed by the combination of the two interactions. That de- 
pends only on the changes in the density matrix p for qubits A 
and B, not on the nature of the controls and the interactions. 
Maps that make the changes in p will be described in the next 
section. 

If there is only the interaction of qubit A with qubit C, the 
concurrence is e~^^*. If there is only interaction of qubit B 
with its control, the concurrence is e~'^^*. If there are both 
and both are made with rotations around the z axis, the con- 
currence is e^'^'^'e^'^^*. If there are both and the interaction 
of qubit B with its control is made with a rotation around the 
X axis, the concurrence is 

C(p±) = imax[0, e-^-**+e"^-**e"^«*+e"^«*-l]. (3.5) 



This concurrence ( 13.51) is zero when 



1. 



(3.6) 



Then the state is separable; it is a mixture of six products of 
pure states: from Eqs.( l3.1l ) and ( I3.3l l 



P'± - ie-i^-*i(l+Si)i(l±5i) 



Vi^-*i(l-I]i)i(]lTSi) 



-TAt-TBt 



1 -r.*l 



1 r A 



2' '2 
■(1+ £3)^(1 -S3) 

■(1 -£3)^(1+53). 



(3.7) 
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The state remains separable at later times; when the sum of 
the exponential decay factors is less than 1, the density ma- 
trix is a mixture in which just a multiple of the density ma- 
trix 1 /4 for the completely mixed state is added to the terms 
of Eq. (13.7b . This change of maximally entangled states to 
separable states can be described without reference to expo- 
nential decay by continuing to use cos (f)A and cos (f)B instead 
of e^'^** and e^^^^. Similar behavior involving exponential 
decay has been observed in more physically interesting and 
mathematically complicated models ISlSlHl- 

IV. MAPS 

The maps that make the changes in the density matrix p for 
qubits A and B could be described in different ways using 
various matrix forms. That is not needed here. Writing p in 
terms of Pauli matrices provides a very simple way to describe 
the maps. For any density matrix 

3 3 3 

^ j = l k=l 3.k=l ' 

(4.1) 

for qubits A and B, the result of the interaction of qubit A 
with qubit C, described by Eq. ( I2.7l i. is that in p, in both the 
and SjSfc terms. 

El — > Si cos(?!)A, E2 — >I]2Cos0a; (4.2) 

the result of the interaction of qubit B with its control is that 
in p 

51 >SiC0S(/)_B, 52 ^S2C0S(/)B (4.3) 

if the interaction is made with a rotation around the z axis; 
and the result of the interaction of qubit B with its control is 
that in p 

52 > S2 cos 0B, S3 > S3 COS (i>B (4.4) 

if the interaction is made with a rotation around the x axis. 
The terms with sin cancel out because there is an equal mix- 
ture of parts with ^ and parts with —0. 

The changes in the state of qubits A and B can be described 
equivalently by maps of mean values that describe the state: 
the result of the interaction of qubit A with qubit C, described 
by Eq. (Ejll, is that 

(El) (El) cos (/.A 

(E2) (E2)C0S(/.A 

(EiSfe) — > (EiSfc)cos</>A 

(EaSfe) (E2Sfc)cos</.A (4.5) 

for fc = 1,2,3; the result of the interaction of qubit B with its 
control is that 

(51) * (Si)c0S(/)B 

(52) > (S2)C0S(/)B 

(EjSi) — > (EjSi)cos0B 

(EjSa) (E,S2)cos0B (4.6) 



for i = 1 , 2, 3 if the interaction is made with a rotation around 
the z axis; and the result of the interaction of qubit B with its 
control is that 

(52) ■> {E2) cos (j)B 

(53) > {E3) COS (j)B 

(EjS2) > (EjS2)cOS0B 

(EjSs) (EjS3)cos0B (4.7) 

for j = 1 , 2, 3 if the interaction is made with a rotation around 
the x axis. 

When (pA and (j>B change over intervals from initial values 
(j>Ai and <f>Bi to final values (j>Af and (j>Bf, the cos0a and 
cos (pB factors in the maps are replaced by cos <j>Af / cos (pAi 
and cos <j>Bf / cos (pBi- If either of these factors is larger than 
1, the map is not completely positive and does not apply to all 
density matrices p for qubits A and B. This happens whenever 
the entanglement increases. It also happens in cases where the 
concurrence (13.21 ) decreases, when one of cos (pA and cos (pB 
increases and the other decreases and there is more decrease 
than increase. The completely positive maps that decrease the 
entanglement have already been described iH]. 

V. RECONCILIATION 

Entanglement being increased by local interactions may 
seem surprising from perspectives framed by experience in 
common situations where it is impossible. Entanglement is 
not increased by a completely positive map of the state of two 
qubits produced by an interaction on one of them. The in- 
teraction will produce a completely positive map if it is with 
a control whose state is initially not correlated with the state 
of the two qubits, as in Eq. ( 12. It . In our examples, that hap- 
pens only when the initial value of is or a multiple of tt. 
Otherwise, the state of the control is correlated with the state 
of the two qubits as in Eq. ( 12.161 ). When a subsystem is ini- 
tially correlated with the rest of a larger system that is being 
changed by unitary Hamiltonian dynamics, the map that de- 
scribes the change of the state of the subsystem generally is 
not completely positive and applies to a limited set of subsys- 
tem states 1 16, 17]. We see this in our examples whenever the 
entanglement increases and in some cases when the entangle- 
ment decreases. 

The map depends on both the dynamics and the initial cor- 
relations. It describes the effect of both in one step. Com- 
pletely positive maps are what you get in the simplest set-up 
where you bring a system and control together in independent 
states and consider the effect of the dynamics that begins then. 
Dynamics over intervals where the maps are not completely 
positive can be expected to play roles in more complex set- 
tings. We should not let expectations for completely positive 
maps prevent us from seeing things that can happen. 

Our perspective is enlarged when we look beyond the map 
and include the dynamics. We can see the dynamics and the 
initial preparation as two related but separate steps. We can 
consider the effect of the dynamics, whatever the preparation 
may be. 
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Local interactions that increase entanglement are com- 
pletely outside a perspective that is limited to pure states. An 
interaction on one of the qubits can not change the entangle- 
ment at all if the state of the two qubits remains pure jlSil . 
The entanglement of a pure state of two qubits depends only 
on the spectrum of the reduced density matrices that describe 
the states of the individual qubits, which is the same for the 
two qubits. If that could be changed by an interaction on one 
of the qubits, there could be a signal faster than light. In our 
examples, the state of the two qubits is pure only when it is 
maximally entangled. In our examples, the spectrum of the 



density matrices for the individual qubits never changes, and 
gives no measure of the entanglement. 
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